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We study the k meson in 2 + 1 flavor QCD witli sufficiently light u/d quarks. Using numerical 
simulation we measure the point-to-point k correlators in the "Asqtad" improved staggered fermion 
formulation. We analyze those correlators using the rooted staggered chiral perturbation theory 
(rSxPT), particular attention is paid to the bubble contribution. After chiral extrapolation, we 
obtain the physical k, mass with 828 ±97 MeV, which is within the recent experimental value 800 ~ 
900MeV. These numerical simulations are carried out with MILC 2-1-1 flavor gauge configurations 
at lattice spacing a ~ 0.12 fm. 
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I. INTRODUCTION 

The so-called k meson (J^ = 0+) is a scalar me- 
son with strangeness. In 2010, Particle Data Group 
(PDG) [l| lists the meson (800) or k with a very broad 
width (550 MeV). A recent analysis [2| gives its mass 
about 75O5Q MeV. Moreover, the resonance of a scalar 
meson is reported 041] to exist in the ttK system with a 
mass of the n meson about 800 MeV. 

Until now, four lattice simulations of the k mass have 
been reported. Prelovsek et al. reported a crude esti- 
mation of the K mass as 1.6 GeV by extrapolating the 
ao mass @. Mathur et al. 0] studied the us scalar me- 
son in the quenched approximation, and obtained the 
value of the k mass to be 1.41 ± 0.12 GeV after remov- 
ing the fitted nrj' ghost. With the dynamical Nf = 2 sea 
quarks and a valence strange quark, the UKQCD Collab- 
oration Q suggested the k mass around 1000—1200 MeV. 
The SCALAR Collaboration 0, [l^ reported full QCD 
simulations on the k meson using the dynamical fermion 
for the light u/d quark and the valence approximation 
for the strange quark, which shows that the / = 1/ 2 
scalar meson has a mass around 1.8 GeV. In Ref. [ll|, 
they performed a quenched QCD calculation using the 
Wilson fermions with the plaquette gauge action, and 
estimated the k mass to be about 1.7 GeV. 

In the presence of 2+1 flavors of Asqtad improved stag- 
gered dynamical sea quarks, we obtained the n mass as 
low as 826±119 MeV in our previous study [l^l, however, 
here we neglected the taste-symmetry breaking and used 
the crude linear extrapolation. It is well known that, in 
the staggered fermion formulation of lattice QCD, due 
to the taste-symmetry breaking there exist many multi- 
hadron states with = 0+ which can proliferate be- 
tween the source and sink of the k correlator. Of spe- 
cial interest for us are the two-pseudoscalar states (i.e., 
bubble contribution) [3, [111 ■ With sufficiently light u 
quark, the k meson propagator is dominated at large 



time distances by these two-meson states. The bubble 
contributions are significantly affected by the unphysical 
approximations which are often used in lattice simula- 
tions [l3, They are expected to disappear in the 
continuum limit. 

In our previous work [l5l - [l8| , we extended the analysis 
of Refs. [1^ 1 examined the scalar meson correlators in 
lattice QCD with the inclusion of the disconnected dia- 
grams, and carried out a quantitative comparison of the 
measured correlators and the predictions from rS^PT. 
Despite the considerable complexity of the scalar meson 
channels with dozens of spectral components, the rSxPT 
provides a strict framework which permits the analysis 
of the scalar meson correlator precisely in terms of only 
a small number of low-energy chiral couplings, which we 
may determine through fits to the data. In this work 
we extend the analyses of Refs. [Tsl - fisj to the k me- 
son, treat the u quark as a valence approximation quark, 
while the valence strange quark mass is fixed to its phys- 
ical value [T9j considering that the ac meson contains a 
strange quark and a light u quark, we perform a series of 
numerical simulations with MILC gauge configurations 
in the presence of 2 -I- 1 flavors of Asqtad improved stag- 
gered dynamical sea quarks, generated by the MILC Col- 
laboration [2^ , and chirally extrapolate the mass of the 
K meson to the physical tt mass using the popular four 
parameter fit with the inclusion of the chiral logarithms. 



II. PSEUDOSCALAR MESON TASTE 
MULTIPLETS 

In Refs. [ll[3, we give a brief review of the rooted 
staggered chiral perturbation theory with particular fo- 
cus on the tree-level pseudoscalar mass spectrum, and 
achieve the rooted version of the theory through the repli- 
cated theory [2l| . 

The tree-level masses of the pseudoscalar mesons 
are [11 [H 
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where x, y are two quark flavor contents which make up, 
b ~ 1, 16 are the taste, fi — m'^/2mq is the low-energy 
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chiral coupling constant of the point scalar current to the 
pseudoscalar field, and the term of a^A^ comes from the 
taste symmetry breaking. The nix and niy arc the two 
valence quark masses in the pseudoscalar meson, and rux 
is the light valence u/d quark mass by convention. 

In this work we investigate degenerate u and d quarks, 
treat u quark as a valence approximation quark, while the 
valence strange quark mass is fixed to its physical value 
nis, thus it will be convenient to introduce the notations 

Mlr^ = AU^ = 2^m^ + a^Ab 

= M.s.b = 2fims + a^Ab (2) 
M](b = ^'^Ki = K'>TT-x + ms) + a^Ab, 

where M(j is the mass of the Goldstonc pion with two 
light valence quark masses, is the mass of the Gold- 
stonc kaon with one valence quark equal to the light va- 
lence quark and one at its physical mass ms, and Ms is 
the mass of a fictitious meson ss in a fiavor nonsinglet 
state (23j with two valence quarks at physical mass m^. 

The isosinglet states (ry and rj') are modified both by 
the taste-singlet anomaly and by the two-trace (quark- 
line hairpin) taste- vector and taste-axial- vector operators 
[23I [2^ . When the anomaly parameter mg is large, we 
obtain the usual result 

</ = + ^MIj, M,,j = 0{ml). (3) 

In the taste-axial-vector sector we have 

M^A = \\MIa + Mh + - Za] 
M^yA = \ [Ml A + Ml A + \^A + Za] (4) 
Zl = {Ml^ - Ml^f - ^(A/J^ - Ml^) + ^S\, 

and likewise for V A, where Sy = a^S'y is the hairpin 
coupling of a pair of taste-vector mesons, and Sa = a^<^A 
is the hairpin coupling of a pair of taste-axial mesons {S'y 
and S 'a a, re tree-level (LO) taste- violating hairpin param- 
eters ^1). 

In the taste-pseudoscalar and taste-tensor sectors, in 
which there is no mixing of the isosinglet states, the 
masses of the rjb and rj'^^ by definition are 

Mlb = Mlb, Mly,b = Mlb- (5) 

In Table m we list the masses of the resulting taste multi- 
plcts in lattice units for our chosen lattice ensemble with 
the taste-breaking parameters Sa and Sy determined in 
Refs. [H [ill. For Goldstone muhiplet (taste P), we 
measured their corresponding correlators and fitted them 
with a single-exponential [1^. Then, using the taste 
splittings in Refs. [l^ we calculated the masses of 
other non- Goldstone taste multiplets. We do not need 
the rj'j masses in this work, hence, we do not list these 
values in Table HI 



TABLE I: The mass spectrum of the pseudoscalar meson 
for MILC medium-coarse (a = 0.12 fm) lattice ensemble with 



/3 = 6.76, 


an^ud = 


0.005, am's 


= 0.05. 






arrix 


taste( B) 


0,71 B 


clKb 


arjB 


o,t1b 




P 


0.1598 


0.3106 


0.1598 


0.4087 


0.005 


A 


0.2342 


0.3546 


0.1831 


0.4332 


T 


0.2688 


0.3784 


0.2688 


0.4624 




V 


0.2971 


0.3990 


0.2832 


0.4755 




I 


0.3198 


0.4162 


0.4434 






P 


0.2233 


0.3291 


0.2233 


0.4087 


0.010 


A 


0.2814 


0.3710 


0.2400 


0.4335 


T 


0.3108 


0.3938 


0.3108 


0.4624 




V 


0.3356 


0.4136 


0.3233 


0.4755 




I 


0.3558 


0.4302 


0.4525 






P 


0.2718 


0.3468 


0.2718 


0.4087 


0.015 


A 


0.3212 


0.3867 


0.2851 


0.4339 


T 


0.3473 


0.4087 


0.3473 


0.4624 




V 


0.3696 


0.4278 


0.3584 


0.4756 




I 


0.3881 


0.4439 


0.4612 






P 


0.3127 


0.3637 


0.3127 


0.4087 


0.020 


A 


0.3565 


0.4020 


0.3235 


0.4345 


T 


0.3802 


0.4231 


0.3802 


0.4624 




V 
V 


U.^UU / 


/i/i 1 




n /i7^7 

U.^ ( ( 




I 


0.4178 


0.4572 


0.4698 






P 


0.3489 


0.3798 


0.3489 


0.4087 


0.025 


A 


0.3886 


0.4166 


0.3572 


0.4356 


T 


0.4105 


0.4370 


0.4105 


0.4624 




V 


0.4295 


0.4550 


0.4195 


0.4759 




I 


0.4455 


0.4701 


0.4782 





Wc should remind the readers that the data in Table U 
are obtained with the valence strange quark mass fixed 
to its physical value. Here and below, we adopt the no- 
tation in Ref. [l9j . the primes on masses indicate that 
they are the dynamical quark masses used in the lattice 
simulations, not the physical masses uiu, md and m^. 



III. THE K CORRELATOR FROM S^PT 

In Refs . [ l5l -[l8j . using the language of the replica 
trick [2^ |25| and through matching the point-to-point 
scalar correlators in chiral low energy effective theory 
and staggered fermion QCD, we rederive the "bubble" 
contribution to the oq channel of Ref. [3], and extend 
the result to the a channel. Here wc further extend these 
results to the k channel. 
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A. Non-singlet k correlator 

To simulate the correct number of quark species, we 
use the fourth-root trick, which automatically performs 
the transition from four tastes to one taste per flavor for 
staggered fermion at all orders. We employ an interpo- 
lation operator with isospin 7=1/2 and = 0+ at the 
source and sink. 



0{x) 



1 



(6) 



where g is the indices of the taste replica, is the num- 
ber of the taste replicas, a is the color indices, and we 
omit the Dirac-Spinor index. The time slice correlator 
C{t) for the k meson can be evaluated by 

X .a.b g,g 

where 0, x are the spatial points of the k state at source, 
sink, respectively. After performing Wick contractions 
of fermion fields, and summing over the taste index [l^ . 
for the light u quark Dirac operator and the s quark 
Dirac operator Mg, we obtain 

= E(-)" (Tr[AC'(a; ,t;Q,0)M-'\x ,t;0,0)]) , 

X 

(7) 

where Tr is the trace over the color index, and x — {x ^t) 
is the lattice position. 

The mass of the k meson can be reliably determined on 
the lattice simulation. However, there exist many mul- 
tihadron states with = 0+ which can propagate be- 
tween the source and the sink. Of special interest for 
us is the intermediate state with two pseudoscalars P\ Pp, 
which we refer to as the bubble contribution (i?) [ij]. 
If the masses of Pi and P2 are small, the bubble term 
gives a considerable contribution to the n correlator, and 
it should be included in the fit of the lattice correlator in 
Eq. ([7]), namely. 



(8) 



where we omit the unimportant contributions from the 
excited states, the oscillating terms corresponding to a 
particle with opposite parity, and other high-order mul- 
tihadron intermediate states. 



B. Coupling of a scalar current to pseudoscalar 

Before we embark on the bubble contribution, here 
we first derive the coupling of a point scalar current 
Sr{x)ur{x) to a pair of the pseudoscalar fields at the 
lowest energy order of the staggered chiral perturbation 
theory (SxPT), where the subscript r in the expression 
Ur{x) is the index of the taste replica for a given quark 
flavor u. The effective scalar current can be determined 



from the dependence of the lattice QCD Lagrangian and 
the staggered chiral Lagrangian on the spurion field M , 
where M. is the staggered quark mass matrix. For n 
Kogut-Susskind (KS) flavors, is a Ann,, x Annr ma- 
trix. 



M-- 



( niuul ' 
rridul ' 



' Ir mudi <■ 
Ir ruddl ( 
Ir msdl ( 



' Ir niusl ( 

Ir mdsl ( 
Ir nissl ( 



^Ir 
'Ir 
Ir 



\ 



\ 



J 



(9) 



where / is a 4 x 4 unit matrix, and Ir is the Ur x unit 
replica matrix. In short, Ai = m(E} I (g) Ir, and 



m = 



^uu ^'ud 

rridu rudd ruds 
rusu rusd niss 



V 



(10) 



is the n X n quark mass matrix. Since the staggered 
chiral lagrangian (CsxPt) is an effective equivalent La- 
grangian for £qcd in low energy limit, the effective cur- 
rent Sr{x)ur{x) is obtained from 



Sr{x)Ur{x) 



dC 



SxPT 



(11) 



where 
CsxPT 



i/^Tr[a^Ea^I]t] - i/i/2Tr[AltS 



is the staggered chiral Lagrangian [2^. We omit the 
high order terms and the terms that are independent 
of A^. /tt is the tree- level pion decay constant [2^ 
is the constant with the dimension of the mass [26 

E = cxp ( ) m. The field $ = ^ -T^ ® 0'' is 

described in terms of the mass eigenstate field cj)'' [2^ . 
where 0'' is a 3 x 3 pseudoscalar matrix with flavor com- 
ponents (ffij J, with flavor /, /', the index of the taste 

replica r, r', and taste b which is given by generators 
= {C5,^^5Ct.,^^^.^.,it.,^I} M- Hence, the field $ 
is AnUr X AnUr pseudoscalar matrix in SxPT |22| . and 
the subscripts u,s denote its valance flavor component. 
Therefore, E is also a Annr x Anur matrix. The Tr is the 
full Annr X Anur trace. Then, the effective current is 129 



~{x)ur{x) — /xTrt[$(x)^ 



(12) 



where the notation Trt stands for the trace over taste, 
and n is the low-energy chiral coupling of the point scalar 
current to the pseudoscalar field $(x). 



C. Bubble contribution to the k correlator 

In this subsection we compute the bubble contribution 
to K correlator in Eq. ([7]) from two intermediate states. 
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From the above discussion, the point scalar current can 
be described in terms of the pseudoscalar field $ by using 
SxPT (Hill, 



Sr{x)Ur{x) = /iTrt [^(x)^]^^^^ , 
Ur{x)Srix) = /iTrt . 



(13) 



For concretencss, the bubble contribution to the k cor- 
relator in a theory with rir tastes per flavor and three 
flavors of KS dynamical sea quarks is [l^] 



r,r — 1 k 
2 "i- f 



bblc 



r,r' — 1 



where the subscript k specifies the bubble contribution 



J 



for K meson. If we consider this identity 
Trt (T'^T'') = Adab, 



we arrive at 



2 16 16 Nf Nf rir Ur 

■■EEEEE E 

a=l b=l 1=1 j=l r,r' = l t,i' = l 



where a,b,a',b' are the taste indices, i,j are the fiavor 
indices, and r,t,r',t' are the indices of the taste replica. 
The Wick contractions result in the products of two prop- 
agators for the pseudoscalar fields. After summing over 
the index of the taste replica, and considering that u, d 
quarks are degenerate, the bubble contribution can be ex- 
pressed in terms of the pseudoscalar propagators {(j)'' (f>^). 



16 



B^^^P^ix) = n^Y^S ('^L(2^)0L(O))con„(^L(a;)'/'L(O))co„n + (^L (2^)0^ (0)) conn (^L (2;)'/>L (0)) c 



b=l 



(0L(^)0L(O))conn(0L(2;)0L(O))c 



+ 



E (</'L(2^)0L(O))conn(0L(x)^L(O)>disc + (0L(a;)0L(O)>conn(0L(-^)0L(O)>di.c 
b=I,V,A I 



2(0L(2:)'^L(O))conn(<^L(a;)<^L(O))di. 



(14) 



The subscript conn and subscript disc stand for the 
connected contribution and the disconnected contribu- 
tion, respectively. The propagators {(jJ^ff)^) for the pseu- 
doscalar field (f)^ for various tastes h are intensively 
studied in Ref. [l^j- The propagators for all tastes 
(/, y, A, T, P) have the connected contributions, while 
only tastes /, V and A have the disconnected contribu- 
tions m. 



('/'L(^)0L(o))Lc = 



-ML) 



(fc2 + Af2J(P + Af2J(fc2 + M2^)' 



UU 



Ml 



3(fc2+M2 )(fc2+Af2 )' 

"3(P + A/|J(P + A/,^,)' 



3fc2 



('^L(2:)'^L(o)) 



disc 



-5y(fc2 + M2j 



(0L(^)'/'L(o))L 



(P + M|J(P + M2J(fc2 + M2^)^ 





■(P+M2J(P+M2^)' 



And likewise for the axial taste (A), we just require 
V ^ A. By plugging in the above mesonic propaga- 
tors (namely, carrying out the Wick contractions) and 
switching to momentum space, Eq. (jl4p can be rewritten 
as 
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1 12 1 1 1 1 

^ (fc + p)2 + M^^ fc2 + ^J2^ + 3 (fc + p)2 + Mf,^ ^2 + ^J2^ " ^ (fc + p)2 + M|.^ fc^ + ^2^ 

<5y fc' + Mjy. fc' + 

(fc + p)2 + M^^. + M|,^, ) (fc2 + Ml ) (F + M2 ^) {k+pr + Ml^ (fc2 + ) (fc^ + Af^^ ) + ^) 

'^y fc' + ^^^v fc' + ^^^A 

{k+pf + Ml^ + M|J(fc2 + A/2J(fc2 + ^2 ^) (fc +p)2 + (fc2 + M2j(p ^ ^,^2J(fc2 + ^^2^) 

g 1 g '^A 1 1 . 

(fc + M|^, + + M2 ^) (fc + (fc2 + M2J(P + M2^) J' ^ 



fe I 6=1 



And likewise for the axial taste {A), we just require 
V A. In Table m we list all these pseudoscalar masses 
needed for the current study. 

In the continuum limit, taste-symmetry is restored, 
namely, 6v = a^5'v ^ 0, (5a = a^5\ 0, then Eq. (fTS|) 
reduces to 

r 16 1 I 

i..(p)^,^E|2..^E(,+,)2 + Mi^P + M^/ 

16 1 r, 1 

Crii- 1 2 1 

[k+pf + Ml^ k^ + M^ ^ (fc+p)2 + M2.^ P + M2^ 

2 1 1 4 1 1 

^3 (/c+p)2 + M|.^ k^ + M^j ~ {k+pY + Mj;^ fc2 + Af|J ■ 

Here the total contribution from pairs of the states with 
mass Mu and Afs is proportional to 

(16n^-4), (16) 

which vanishes when = 1/4. The negative thresh- 
old has nicely canceled out the unphysical threshold KS. 
The surviving thresholds are the physical -qK. 

IV. SIMULATIONS AND RESULTS 

We use the MILC lattices with 2-1-1 dynamical flavors 
of the Asqtad-improved staggered dynamical fermions, 
the detailed description of the simulation parameters can 
be found in Refs. [lO, 0, 1131 . We analyzed the k correla- 
tors on the 0.12 fm MILC ensemble of 520 24^ x 64 gauge 
configurations with bare quark masses am'^^ = 0.005 
and am\ = 0.05 and bare gauge coupling lO/^^ = 6.76 
the lattice spacing oT^ = 1. 679^^3 GeV, which has a 
physical volume approximately 2.5 fm. The mass of the 
dynamical strang e q uark is close to its physical value 
ams = 0.0344 [23, [231. The masses of the u and d quarks 
are degenerate. In Table HI we list all the pseudoscalar 
masses used in our fits with the exception of the masses 



MriA ; ^^T]'^ 1 MrjY , and M,,^ . Those masses can vary with 
the fit parameters Sa and 6v- 

For the light u quark Dirac operator M^ and the s 
quark Dirac operator Ms , we measure the point-to-point 
quark-line connected correlator which is described by 
Eq. ([7]). We use the conjugate gradient method (CG) to 
obtain the required matrix element of the inverse fermion 
matrix M„ and Ms- 

In order to improve the statistics, we place the source 
on all the time slices tg = 0, • • • ,T — 1, therefore, we 
perform T = 64 inversions for each configuration and av- 
erage these correlators. Note that the time extent of our 
lattices is more than twice the spatial extent. The rather 
large effort to generate propagators allows us to evaluate 
the correlators with high precision, which is important 
to extract the desired k masses reliably. 

Since the k meson contains a strange s quark and a 
light u quark, we should treat the u quark as a valence 
approximation quark, while the valence strange quark 
mass is fixed to its physical value [l^ . The physical value 
of the strange quark mass of the lattice ensemble used in 
the present work has been precisely determined by MILC 
simulations [20| . namely, arris = 0.0344, where a is the 
lattice spacing. 

The propagators of the k meson are calculated with the 
same configurations using five u valence quarks, namely, 
we choose am^ = 0.005, 0.01, 0.015, 0.02 and 0.025, 
where rrix is the light valence u quark mass. In order 
to obtain the physical mass of the k meson, we then per- 
form extrapolation to the chiral (physical tt mass, ob- 
tained from PDG) limit guided by chiral perturbation 
theory. The correlators of the tt, K meson and fictitious 
meson ss are also measured with the same configurations 
for calculating the pseudoscalar masses in Table HI 

Figure [1] shows the k propagators with five different 
light valence u quark masses and their predicted bubble 
contributions. For our chosen MILC configurations used 
in the present study we obtain the positive predicted bub- 
ble contributions, since the small negative contribution 
in the bubble term is most likely outweighed by the pos- 
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FIG. 2: The k masses as a function of the minimum time 
distance. The effective mass plots will be a plateau in time 
range 9 < Anin < 11- 



FIG. f : The k propagators for five valence it quarks. Overlaid 
on the data are their corresponding predicted bubble contri- 
butions. 

itivc contribution as is discussed for the og correlators in 
Ref. [l^. We add a constant 7.8e — 6 to all data points 
and the corresponding bubble terms in y axis just for 
good visualization. Figure [T] clearly shows that the pre- 
dicted bubble contributions dominate the k propagators 
after t > 12. 

For staggered quarks the meson propagators have the 
generic single-particle form. 

C{t) ^ ^«e'™'* + A{-iye-<* + it^Nt- t), 

i i 

(17) 

where the oscillating terms correspond to a particle with 
opposite parity. For the k meson correlator, we consider 
only one mass with each parity in the fits of Eq. (1171) . 
namely, in our concrete calculation, our operator is the 
state with spin-taste assignment I ® I and its oscillating 
term with spin-taste assignment 7075 ® 7075 [i3 • From 
the aforementioned discussion, we must consider the bub- 
ble contribution. Therefore, all five k correlators were 
then fit to the following physical model 

a(t)=C~(t) + B„(t), (18) 

here 

where the and 6„ are two overlap factors, and the 
bubble term in the fitting function Eq (jl8p is given in 
momentum space by Eq. (jisp . The time- Fourier trans- 
form of it yields 3^,(1:,). 

This fitting model contains the explicit n pole, together 
with the corresponding negative-parity state Ka and the 
bubble contribution. There are four fit parameters (i.e., 



, Mka I ^Ka I ^i^d 6k ) for each k correlator with a given 
valence u quark mass mx ■ The bubble term B^^ (t) was pa- 
rameterized by three low-energy couplings /i, Sa, and Sy- 
In our concrete fit, they were fixed to the values of the 
previous MILC determinations [l^. The taste multiplet 
masses in the bubble terms were fixed as listed in Table ID 
The sum over intermediate momenta was cut off when 
the total energy of the two-body state exceeded 2.0/a or 
any momentum component exceeded 7r/(4a). We deter- 
mined that such a cutoff gave an acceptable accuracy for 
t > 8. The lightest intermediate state in bubble term is 
ttK. Therefore, this fitting model can remove a number 
of unwanted ttK states with different tastes and slightly 
different energies. 

For arrix = 0.005, the effective mass plots of the k 
meson are shown in Figure [2l We find that the effec- 
tive K mass suffers from large errors, especially in larger 
minimum time distance regions. To avoid possible large 
errors coming from the data at large minimum time dis- 
tance Dmin, we fit the effective mass of the k meson only 
in the time range 9 < i^min < H, where the effective 
masses are almost constant with small errors. 

In our fit, five k propagators were fit using a minimum 
time distance of 10a. At this distance, the contamination 
from the excited states is comparable to the statistical 
errors, we can neglect the systematic effect due to excited 
states, therefore we can extract the mass of the k meson 
efficiently. 

The fitted masses of the k correlators are summarized 
in Table |lll The second block shows the masses of the 
K meson in lattice units, and Column Four shows the 
time range for the chosen fit. As a consistency check, we 
also list the fitted masses of their corresponding negative 
parity state Ka in Column Three. We can note that the 
fitted values of the pseudoscalar meson Ka masses are 
consistent with our calculated values in Table |I] within 
small errors. Column Five shows the number of degrees 
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TABLE II: The summary of the results for the fitted k 
masses. The second block shows the k masses in lattice units. 
The third block shows the fitted Ka masses. 



arrix 




aMKA 


Range 


xVdof 


0.005 


0.545(26) 


0.3568(25) 


10 


-25 


12.7/12 


0.010 


0.590(23) 


0.3719(17) 


10 


-25 


11.8/12 


0.015 


0.625(21) 


0.3868(14) 


10 


-25 


11.2/12 


0.020 


0.653(19) 


0.4031(12) 


10 


-25 


11.0/12 


0.025 


0.673(18) 


0.4182(11) 


10 


-25 


11.4/12 



of freedom (dof) for the fit. 

In order to obtain the physical mass of the k meson, 
we carry out the chiral extrapolation of the k mass to 
the physical tt mass using the popular three parameter 
fit with the inclusion of the next-to-next-to-leading order 
(NNLO) chiral logarithms. The general structure of the 
pion mass dependence of to^ can be written down as 

"^K = Co + C2m^ + cam^ -I- 0477^ In(m^), (19) 

where cp, C2, C3 and C4 are the fitting parameters, and the 
fourth term is the NNLO chiral logarithms. 

We obtain the physical tt mass from PDG [l| and use it 
as a chiral limit. In Figure[3l we show how physical value 
TOk is extracted, which gives x^/dof = 0.16/1. The blue 
dashed line in Figure [3] is the linear extrapolation of the 
mass of the k, meson to the physical pion mass ttItt- The 
chirally extrapolated k mass = (828±97) MeV, which 
is in good accordance with the result in our previous 
study on a MILC medium coarse lattice ensemble [l^ . 
The cyan diamond in Figure [3] indicates this value. 

Using the fitting model in Ref. [1^, we extract kaon 
masses. In Figure |31 we display to^, fn-n, and ttIt^j^k 
in lattice units as a function of the pion mass m-^. We 
observe that, as the valence quark mass increases, ttK 
threshold grows faster than k mass and, as a consequence, 
■nK threshold is higher than k mass for large pion mass 
(about arrix > 0.151), and ttK threshold is lower than k 
meson for small pion mass (about anix < 0.151). Hence, 
the K meson can decay on our lattice for small quark 
mass, but for large quark mass, the decay k — ttK is no 
longer allowed kinematically, which is in good agreement 
with the results in Ref. (28| . 

To understand the effects of these bubble contributions 
(or taste breakings), in the present study we also fitted 
our measured kappa correlators without bubble terms. 
The fitted results are tabulated in Table IIIIl From Ta- 
ble In] and Table IIIIl we can clearly see that the bubble 
terms contribute about 2% — 5% differences for the k 
masses. 



V. SUMMARY AND OUTLOOK 

In the present study we have extended the analyses of 
the scalar mesons in Refs. and derived the two- 
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FIG. 3: Characteristics of 7TiK,mK,?Ti^ and m,r+K in lattice 
units as a function of the pion mass. The chiral limit is ob- 
tained at the physical pion mass rrin- 



TABLE III: The summary of the results for the fitted k 
masses without bubble contributions. The second block shows 
the K masses in lattice units. The third block shows the fitted 
Ka masses. 



arrix 




uMka 


Ran 


ge 


xVdof 


0.005 


0.572(19) 


0.3565(25) 


10- 


25 


12.1/12 


0.010 


0.615(19) 


0.3717(17) 


10- 


25 


11.0/12 


0.015 


0.646(19) 


0.3876(14) 


10- 


25 


10.4/12 


0.020 


0.669(18) 


0.4030(12) 


10- 


25 


10.2/12 


0.025 


0.687(18) 


0.4181(11) 


10- 


25 


10.7/12 



pseudoscalar-meson "bubble" contribution to the k cor- 
relator in the lowest order S^PT. We used this physical 
model to fit the lattice simulation data of the point-to- 
point scalar k correlators for the MILC coarse (a « 0.12 
fm) lattice ensembles in the presence of 2 -I- 1 flavors of 
Asqtad improved staggered dynamical sea quarks, gener- 
ated by the MILC Collaboration [Uliil. We treated the 
light u quark as a valence approximation quark, while 
the strange valence s quark mass is fixed to its physical 
value, and chirally extrapolated the mass of the k me- 
son to the physical pion mass. We achieved the physical 
mass of K meson with 828 ± 97 MeV, which is very close 
to the recent experimental value 800 ~ 900 MeV. Prob- 
ably, it may be identified with the k meson observed in 
experiments. 

Most of all, we note that k meson is heavier than ttK 
threshold for enough small u quark mass. Therefore, it 
can decay on our lattice for small quark mass. This pre- 
liminary lattice simulation will stimulate people to study 
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the decay mode k — )■ ttK. We are beginning lattice 
study of this decay channel with isospin representation 
of / = 1/2. 

Since the appearance of the bubble contribution is a 
consequence of the fermion determinant, an analysis of 
the K correlator in this work also provides a direct useful 
test of the fourth-root recipe. The bubble term in S^PT 
provides a useful explanation of the lattice artifacts in- 
duced by the fourth- root approximation [l6j . The ar- 
tifacts include the thresholds at unphysical energies and 
the thresholds with negative weights. These contribu- 
tions are clearly present in the k channel in our QCD sim- 
ulation with the Asqtad action at a ~ 0.12 fm. We find 
that the "bubble" term must be included in a successful 
spectral analysis of the k, correlator. The rSxPT predicts 
further that these lattice artifacts vanish in the contin- 
uum limit, leaving only physical two-body thresholds. It 
would be nice to be able to investigate whether this ex- 
pectation is ruled out in lattice simulations at smaller 
lattice spacing. 

In this work we reported our preliminary results on 
one lattice ensemble, more physical one should be in the 
continuum limit. We are beginning a series of numerical 
simulations with the MILC fine and super-fine lattice en- 
sembles. Furthermore, the kappa meson is a resonance, 
i.e. a state with a considerable width under strong in- 



teractions. In order to map out "avoided level crossings" 
between the resonance and its dec ay p roducts in a finite 
box volume proposed by Liischer [29| . we are beginning 
to measure the itK scattering K+t: — > K-'rn channel and 
the K i^-fTT, K+t: — ^ k cross-correlators in addition to 
the K — > K correlator [30l - [33l | . Hopefully, we can reliably 
obtain the resonance parameters of the k resonance. 
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